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Abstract 

The unification of gauge coupling constants in the minimal supersymmetric model 
(MSSM) is unaffected at the one-loop level by the inclusion of additional mass- 
degenerate SU(5) multiplets. Perturbativity puts an upper limit on the number of 
additional fields. We analyse the evolution of the gauge coupling constants in all 
models satisfying these criteria using two-loop (3 functions and including low energy 
threshold effects. We find that similarly to the minimal supersymmetric model 
(MSSM) unification takes place within the theoretical and experimental errors. The 
dominant proton decay mode is more suppressed in all extended models as opposed 
to the MSSM due to renormalization group effects. However, the prediction for the 
bottom to r mass ratio becomes worse in all models under consideration. 





The prediction of the strong coupling constant, a s , within the framework of the minimal 
supersymmetric model (MSSM) assuming gauge coupling unification without any inter¬ 
mediate scale is in acceptable agreement with experiment Q. In addition, the unification 
of r and bottom Yukawa couplings is quite promising 

Despite these successes of supersymmetric grand unified theories (SUSY GUTs) based on 
SU(5)|| there remain still some problems. For example, the prediction for the down type 
quarks from Yukawa unification of the first two generations is off by a factor / (■ m e m s ) 

0(10). Maybe the most sever challenge is the so-called doublet/triplet problem^ of giving 
the colored Higgs triplet a mass of the order of Mq ut while retaining the Higgs doublets 
responsible for the electro-weak symmetry breaking at the electroweak scale. There have 
been many attempt to try and solve these problems all of which have one thing in com¬ 
mon: they require the introduction of new fields. Thus, an extension of the particle 
content of the MSSM is inevitable and one might ask whether some of these new particles 
are present in the low energy effective theory. This idea has already been explored in 
ref. H |j] and will be reconsidered here in view of the improved experimental limits on 
the gauge coupling constants. In addition, we will relax the constraint for perturbativity 
and as a result obtain models that have not been considered before. 

From the severity of the doublet/triplet problem we know how hard it is to construct a 
model with a large mass hierarchy between different SU(3)®SU(2)®U(1) members of the 
same SU(5) multiplet. Therefore, we can expect that the most natural extensions of the 
MSSM will come as complete SU(5) multiplets in the low energy effective theory. It is a 
simple exercise to show that all such extensions will preserve gauge coupling unification 
at the one-loop level. 

The renormalization group equations (RGEs) for these models can be written as 

n ry ■ 

— G- = of [A + 0ijf3ij + 0(a 2 )} . (1) 

Here, t = (27r) _1 In (scale), the indices i,j = 1,2,3 refer to the U(l), SU(2) and SU(3) 
gauge group and summation over twice occurring indices is assumed. Furthermore, the 
one-loop (3 functions for the gauge couplings are[[fj] 
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where the three contributions to (3i come from the gauge sector, the Higgs doublets (in the 
MSSM the number of Higgs doublets, N H = 2) and the contribution of complete SU(5) 
multiplets. The contributions of the gauge/gaugino sector to /?* are non-universal since 
some gauge bosons and their superpartners acquire a mass via the Higgs mechanism while 
others stay massless due to gauge invariance. The contributions of the Higgs bosons are 
also non-universal because the doublets are responsible for the electro-weak symmetry 
breaking and should have a mass of the order of 100 GeV while the Higgs triplets have 
to acquire a mass at Mqut in order to sufficiently suppress the rate for proton-decay. 
This doublet/triplet splitting^ is solely motivated by experiment and has no satisfying 
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theoretical solution yet. The last term which also contains the contribution of Nq = 3 
generations of quarks and leptons is universal for all three couplings. (Note, that one 
family of quarks and leptons can be embedded in a 5 = d c { 3,1,2/3) © Z(l, 2, —1) and a 
10 = q( 3, 2,1/3)©« C (3,1, —4/3)©e c (l, 1, —2); the numbers in brackets indicate the trans¬ 
formation properties under the SU(3), SU(2) and U(l) gauge symmetries, respectively.) 
The reason is that the inclusion of a full SU(5) multiplets with a mass, m, does not break 
the SU(5) gauge symmetry and should yield a universal contribution to all three (5 func¬ 
tions at any scale above m at the one-loop level. The contributions of the extensions of 
the MSSM can be written as 


i3 x = '£n$), ( 3 ) 

$ 

where the sum is over all SU(5) multiplets $. The values of T(<f>) = d(<h)C' 2 ( < l ) )/r are 
listed in Table [I] for the four smallest representations of SU(5). Here, r = 24 is the number 
of generators of SU(5) and C2(4>) [d( < h)] is the quadratic Casimir operator [dimension] of 
the SU(5) representation $||. By imposing gauge coupling unification at Mqut he. 

«GUT = OCi(Mqu T ) = (^(Mtut) = oi 3 (Mq ut) , (4) 

and solving eq. [l] to first order in perturbation theory we obtain 

1 M gut _ a^{m z ) - a^irriz) 

h) = tt” 111 -—- a - a -— °- 6 > 

27t m z /?i - p 2 

«3 1 (m z ) = Qta 1 (m z ) + t 0 ((3 3 - # 2 ) ~ 8 , (5) 


where we have used 

5 ^em 
3 cos 2 9 m g ’ 
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( 6 ) 


and a~^ = 127.9 and sin 2 9 = 0.23190 as the low energy input values. The prediction 
of the strong coupling constant, a s (m z ) = a 3 (?n z ) ~ 0.125 is in quite good agreement 
with the world average a s (m z ) = 0.117 ± 0.005 0. In deriving eq. |] we have assumed 
that there is no intermediate scale but it also holds in the case of a widely spread particle 
spectrum as long as the members of the different SU(5) multiplets lie close together. 

Note, that the right hand side of eq. [5] is independent of /3 X . This means that any 
extension of the MSSM by full SU(5) multiplets will maintain the property of gauge 
coupling unification at one loop. The unification scale Mgut — 2 x 10 16 GeV remains 
also unchanged and for the unified gauge coupling we obtain 


aT 1 ~ 24 - t 0 (5 x . 


(7) 


By requiring that the right hand side of eq. |7] is larger than zero we find (3 X ^ 4.5 but 
maybe models with /3 X = 5 are still acceptable due to higher order corrections or threshold 
corrections etc. and shall be included into our considerations. 
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Table 1: The contribution of various SU(5) multiplets <f> to (3 X . 

In order to derive a viable model, we have to impose additional constraints. The can¬ 
cellation of triangle anomalies implies that complex representations only occur in pairs. 
Thus, there are four types of extensions satisfying the above requirements 

• n additional pairs of 5 and 5, where n = 1,2, 3,4, 5, 

• one additional pair of 10 and 10, 

• n additional pairs of 5 and 10 where n — 1,2, 

• one additional adjoint representation, 24 = g( 8,1, 0)©w(l, 3, 0)©6(1,1, 0)©x(3, 2, 5/3)© 
x(3, 2, —5/3). 

Experimental lower limits on the additional particle masses can be satisfied by adding 
explicit clirac or majorana mass terms into the superpotential such as 

W D /i 5 55, /TiolOlO, ^ 24 24 2 , (8) 

allowed by gauge invariance in models of type 1, 2 and 4. Models of type 3 correspond 
to the MSSM with four or five generations. For the invisible width of the Z boson at 
LEP experiments we know that the number of (almost) massless neutrinos is 3 and a 
mechanism has to be introduced in order to give mass to the additional neutrinos larger 
than about m z /2. Also, a lower limit on the mass of an additional lepton of m T > > m z /2 at 
LEP and a lower limit on the mass of an additional down type quark of mv > 85 GeV at 
CDF have been established|9j]. This implies that the Yukawa couplings for the additional 
fermions are bound from below since no explicit gauge invariant mass term exists. On the 
other hand, there is an upper limit on the masses from the infra-red fixed-point behaviour 
of the Yukawa couplings leaving only a very constrained region in parameter space. The 
four generation model has been studied recently and found to be quite constrained but 
could still be feasible if a right-handed neutrino is introduced to raise the mass of the 
left-handed neutrino above the experimental bounds|l(J. The five generation model is 
even more constrained and might already be ruled out by present data. However, the 
model with additional 5, 5,10,10 with the possibility for a explicit dirac mass for all 
the additional particles is still allowed. It also does not require any additional fields to 
generate mass for the unseen neutrinos. 

The assumption that all the members, 0, of one SU(5) multiplet, $, are mass degenerate 
is protected by gauge invariance. It acquires corrections below M GUT where the gauge 
symmetry is broken through one-loop RG evolution 

> ( 9 ) 
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Tabic 2: The (5 functions for the additional fields, 0 


where the (3^ are listed in Table. |2j. The solution of eq. |9] can be written as 

/v(/v) = n f q gut \ 

dGVT J 


( 10 ) 


where we have assumed that all possible additional couplings of the superpotential are 
much smaller than «gut and can be neglected in the evolution of the ads. Furthermore, 
we assume that the complex (real) multiplets are odd (even) under R-parity so that proton 
decay mediating interactions are forbidden. We see that the splitting between the various 
members of an SU(5) multiplet is quite considerable and gives rise to significant threshold 
corrections. If we decouple the fields 0 from the RGEs at we obtain the improved 
one-loop formulas 


cq 1 (m z ) = a G u T + t 0 A - A 


MSSM 

i 


--^-J2 N 4> ln —A,* , 

2ti m 7 


( 11 ) 


where A,0 is the contribution of the field, 0 to A listed in Table 


11]. Furthermore, 
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(12) 


is the number of the fields 0. The MSSM threshold corrections, A^ SSM , studied in ref. [[12] 
raise the predicted value of a s (m z ) by about 10% and spoil the success of the GUT 
prediction to some degree. 

For the computation of A^ SSM we have assumed that the effects of SUSY breaking are 
parameterized in the standard fashion by including explicit soft SUSY breaking terms 
assumed to be universal at M GUT . Thus, at M GU t we have only four independent soft 
SUSY parameters: the coefficients multiplying the trilinear and quadratic terms of the 
superpotential, A and R, which are irrelevant for the mass spectrum and, thus, for the 
value of the gauge couplings at the leading log level, the mass for all spin 0 particles mo, 
and the gaugino mass parameter, mm. (We refrain from a more sophisticated treatment 
of the SUSY threshold corrections since here we only care about the changes of extended 
models with respect to the MSSM.) From the last two parameters we can derive the full 
low energy mass spectrum via RG evolution. 


In deriving eq. [TT] we have also made the assumption that the mass splitting between spin 
0 and spin 1/2 components due to soft SUSY breaking mass terms can be neglected. This 
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Table 3: The one-loop 0 functions for the SU(3)<S>SU(2)<g>U(l) gauge couplings due to 
the individual representations 0. Note that by summing over a complete SU(5) multiplet 
we recover the universal coefficients of Table 1 


is motivated by the observation that the Higgs mass parameter // has to be of the order 
of the top squark mass for correct radiative electro-weak symmetry breaking. However, 
from eq. |l0| we see that //^(/x^) /x 2 (m z ) if 0 carries color. Also, /x 2 (m z ) is additionally 

suppressed by the effect of the large top Yukawa coupling, a t . We will return to this point 
in the context of radiative electro-weak symmetry breaking. 

Additional effects on the unification of the gauge couplings arise because the universality 
of 0 A ' in eq. 0 is violated at the two-loop level. The two-loop 0 functions of these models 


can easily be derived from ref. Jl3| 
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Here we have not separated out the contributions of the individual representations since 
the two-loop threshold corrections analogous to the one-loop corrections of eq. 
negligible. 
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are 


In minimal SU(5) SUSY-GUTs the down-type quark fields d and the left-handed lepton 
fields l are embedded in one representation and as a result the r and bottom Yukawa 
couplings are unified at Mqut- In the limit of negligible r and bottom Yukawa couplings 
( i.e. a T , a b <C ckgut) we can write in the one-loop approximation 


m b (m z 


m-r m 7 


= n,: 


(ag)) 


«gut 


a. 


(#-#)/& 


exp 


1 

2 




a t dt 


(14) 


where /3 l b = (7/18,3/2,8/3) and 0* = (3/2, 3/2,0). It is easy to see that m,b/m T increases 
with «gut and hence also with /3 X . This increase can be compensated by an increase 
in the top Yukawa coupling, a t . However, a t quickly approaches its IR fixed-point a t — 
(8/9)ct s (near Mqut: a t — (44/27)q:gut)- Hence, any significant increase of the value 
of the integral in eq. [b| can only come from the integration close to Mqut and requires 
unperturbatively large values of a t (M G ut)- Thus, in our numerical work we have used 
«0 Mgut) = c^gut in order to obtain a natural upper limit for the integral in eq. 14 
and, hence, a natural lower bound on mb/rn T for a particular model assuming r-bottom 
Yukawa unification. 

Another important constraint on SUSY-GUT models comes from non-observation of 
proton-decay. In the MSSM the dominant decay proceeds via p — > K + z014] by dressing 
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Table 4: The one-loop (left column) and two-loop (right column) results for the ratio, 
r = Mq UT /Mqu f M ; Aa s /a s [in %], g G ut, the ratio mb/m T . and the ratio r A ’/T MSSM . We 
have chosen m 0 (M G ut) = 200 GeV, /^(Mgut) = 1 TeV for all additional fields, </>, and 
cp(Mgut) = «gut- Furthermore, we have fixed mi / 2 such that Mg(Mg) = 200 GeV. 

of the dimension 5 operators obtained from the non-renormalizable term of the superpo¬ 
tential 


W NR = T(qq)(ql) , 


(15) 


where the SU(2) indices are contracted inside the brackets and flavor and color indices 
are omitted. At the one-loop level we find 


T oc 


M, 




GUT 


' aiK ) 1 

«GUT 


(PL+ft-bt-biVPi 


(16) 


ffere, we have defined f3 l d = f3 l b , (3 l u = (13/18,3/2,8/3), 7 * = (1/10,3/2,8/3), and = 
(1/2, 3/2,0). Of course, a precise determination of the value of T is problematic without 
a selfconsistent model for the origin of the Yukawa couplings for all three generations, 
ffowever, we do not worry about any overall factors since we are only interested in the 
ratio r x /T MSSM . 

In Table [l| we have summarized our results for the MSSM and 8 extended models char¬ 
acterized by IV 5 , N\ 0 , and i\T 24 . We have chosen the somewhat large values of = 1 TeV 
for all additional fields (f) in order to exploit as many models as possible. The parame¬ 
ter independent lower limit on the gluino mass of Mg > 100 GeV has been established 
from direct particle search at CDF [[15). However, stronger limits can be derived from the 
chargino/neutralino search by imposing GUT constraints and we chose Mg = 200 GeV in 
order to safely avoid all the present bounds. 

The first row corresponds to the MSSM where we have = 0.124 and M/fufU = 

2.3 x 10 16 GeV for our choice of parameters. In the different columns we present the 
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result for the unification scale divided by the two-loop MSSM value and denoted by r, 
the relative change in the prediction of a s with respect to the two-loop MSSM value, 
the unified gauge coupling, «gut, and the ratio of rrib(m z ) to m T (m z ). The two values 
correspond to the results obtained by using one-loop and two-loop (3 functions. We see 
that already for the MSSM the value of nrib/m T for a t (M GUT ) = ckgut is slightly above its 
experimental value of rrib/m T ~ 1.6 but can still be brought in agreement with experiment 
by a modest increase of a t or by choosing = 0(a t )|2|. The situation becomes more 
problematic in all extended models where this ratio increases even more. 


Furthermore, we see that the models 4, 7, 9 and 10 become non-perturbative at Mgut- 
This scenario of non-perturbative unification was already advocated in ref. |16| in non- 


SUSY models and in ref. [[[7] extended to SUSY models as being particularly attractive. 
The reason is that the dependence of cq on a G UT in eq. [□] vanishes in the large «gut 
limit. In the presence of explicit mass terms this argument becomes more complicated. 


From eq. [lb| we see that the low energy particle spectrum becomes non-predictive when 
ckgut vanishes. However, we will still be able to use perturbation theory to predict certain 
trends in these models. For example, at the one-loop level the gluino mass satisfies 


. . a s (M 5 ) 

M- g {M - g ) = — —mi/2 


«GUT 


(17) 


and thus we find M g (M g ) -C m \/2 in all models with non-perturbative gauge unification. 
On the other hand, we have for the squark masses 


M l = Cm l/2 + ..., (18) 

where C ~ 6 — 0.8 /3 X 2 and we have dropped irrelevant terms. Thus, by imposing 

M g (M g ) = 200 GeV we find typically that all the squark masses have to be in the TeV 
region or higher depending on «gut- Such a large hierarchy between the squark masses 
and the electro-weak scale requires fine-tuning and is problematic as we will see in the 
following. In the limit of small gluino mass the soft SUSY breaking top squark mass 
parameters, M~ ( P = L,R) and the Higgs mass parameter, M H2 , are closely related by 
the solution to their RG equations 


£(ra z ) 

S(Mgut) 



(19) 


where we have defined £ = + M~ + Mjj 2 . In the limit of large a t the right hand side 


of eq. [19] vanishes and we find 


M 2 H2 


~ Mf + Mr 

tL tR 


> 


(1 TeV)' 


( 20 ) 


The correct clectro-weak symmetry breaking then implies that we have to fine-tune /x and 
M 2 H2 over several orders of magnitude in order to satisfy 


J-4 + y = 0 (ml). 


( 21 ) 


This heavy SUSY particle spectrum and the associated fine-tuning is a generic problem 
in all models with non-perturbative unification over the entire parameter space. 
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Tabic 5: The prediction of the low energy gauge coupling constants for various boundary 
conditions at M GUT in model 10 


It is interesting to note that model 10 automatically provides a solution to this problem. 
Here, the gluino can acquire a mass via mixing with the octet of the 24 dimensional 
representation || and as a result the mass parameter Mg(Mg) can be much smaller than 
the lower bound on the gluino mass. Such a dirac mass term can be compatible with 
R symmetry ||18|| and was already advocated in ref. |l9j as the origin of the gluino mass. 
Furthermore, a model with an approximate R symmetry naturally yields large value of 
Higgs vacuum expectation values, tan f3, which explains the large ratio of rn t to mb- In 
the MSSM this can only be achieved by fine-tuning [I 


In order to see how far we can trust perturbation theory in the case of non-perturbative 
unification we will modify the boundary condition of eq. [4] for model 10. In Table |5| we have 
listed the prediction for a em (m z ), sin 2 g, and a s (m z ). We have used M G ut = 1-6 x 10 16 
and o G ut = 3.2. 


We find indeed that if Mqvt is close to a Landau pole any potentially large threshold 
corrections get washed out via RG evolution into the perturbative region below M G ut- 
We see that changes of the GUT input parameters by a factor of 10 only changes the low 
energy values by typically a few % (the exception is the case where a 3 
loti, a 2 because here, the splitting of the 24 due to eq. 0 is only mild). As a result, a theory 
that is non-perturbative at M GG t might still allow for reliable low energy predictions. 
In fact, in model 10 the three gauge couplings are completely fixed by only one input 
parameter, M G ut, rather than two as in the MSSM. However, it is somewhat surprising 
that the predictability of the model is not so much limited by non-perturbative effects 
at M G ut but rather by the low energy threshold corrections in eq. 0. The reason is 
that the particle spectrum derived in eq. 0 becomes non-predictive in the large « G ut 
limit and can only provide us with a qualitative understanding. In general, we expect 
colored particles to be much heavier than the electro-weak scale and to decouple at a 
higher scale. This predicts an decrease of a. 3 l due to eq. 0 which for our particular 
choice of parameters and using our two-loop RG approach results in an overestimate of 
a s by about 30%. This number should not be interpreted as a prediction but rather as an 
indication for how the unpredictability of the particle spectrum feeds into the prediction 
of the coupling constants. 



















In summary, we have investigated 10 extensions of the MSSM with additional SU(5) 
multiplets. We find that in all models where ckgut < 1 the prediction of a s only changes 
by a few %. Furthermore, the ratio mb/m T increase and the proton decay rate decreases 
with increasing «gut- We point out that models with non-perturbative unification predict 
in general a large squark to gluino mass ratio which leads to a fine-tuning problem once 
we impose the experimental lower limit on the gluino mass. The exception is the model 
with an additional 24 dimensional representation where the gluino can have a dirac mass. 
However, SUSY threshold corrections to a s are positive and large due to large multiplet 
splitting and in general raise a s significantly above the experimental limit. 
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